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ABSTRACT 



The evolution of intermediate mass stars at very low metallicity during their 
final thermal pulse asymptotic giant branch phase is studied in detail. As rep- 
resentative examples models with initial masses of 4 M and 5 M with a metal- 
licity of Z = 0.0001 ([Fe/H] ~ —2.3) are discussed. The ID stellar structure and 
evolution model includes time- and depth dependent overshooting motivated by 
hydrodynamical simulations, as well as a full nuclear network and time-dependent 
mixing. Particular attention is given to high time and space resolution to avoid 
numerical artefacts related to third dredge-up and hot-bottom burning predic- 
tions. The model calculations predict very efficient third dredge-up which mixes 



the envelope with the entire intershell layer or a large fraction thereof, and in 
some cases penetrates into the C/O core below the He-shell. In all cases primary 
oxygen is mixed into the envelope. The models predict efficient envelope burning 
during the interpulse phase. Depending on the envelope burning temperature, 
oxygen is destroyed to varying degrees. The combined effect of dredge-up and 
envelope burning does not lead to any significant oxygen depletion in any of the 
cases considered in this study. The large dredge-up efficiency in our model is 
closely related to the particular properties of the H-shell during the dredge-up 
phase in low-metallicity very metal poor stars, which is followed here over many 
thermal pulses. During the dredge-up phase, the temperature just below the con- 
vective boundary is large enough for protons to burn vigorously when they are 
brought into the C-rich environment below the convection boundary by the time- 
and depth dependent overshooting. H-burning luminosities of 10 5 to ~ 2 x 10 6 L 
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are generated. C and, to lesser degree 0, is transformed into N in this dredge- 
up overshooting layer and enters the envelope. The global effect on the CNO 
abundance is similar to that of hot bottom burning. If the overshoot efficiency 
is larger, then dredge-up H-burning causes a further increase in the dredge-up 
efficiency. After some thermal pulses the dredge-up proceeds through the He- 
shell and into the CO core beneath. Then neutrons may not be released from 
13 C in radiative conditions during the interpulse phase because of the scarcity 
of a-particles for the 13 C(a,n) 16 reactions. Conditions for the s-process are 
discussed qualitatively. The abundance evolution of H, He, C, N, O and Na is 
described. Finally, the model predictions for sodium and oxygen are compared 
with observed abundances. The notion that massive AGB stars are the origin of 
the O-Na abundance anti-correlation in globular cluster giants is not consistent 
with the model predictions of this study. The abundance of the C-rich extremely 
metal poor binaries LP 625-44, CS 29497-030 and HE 0024-2523 is discussed. 

Subject headings: stars: AGB and post- AGB — stars: evolution — stars: abun- 
dances — stars: Population II — nuclear reactions, nucleosynthesis, abundances 

1. Introduction 

During their final evolutionary phase intermediate mass stars (IMS) proceed through 
the thermal pulse Asymptotic Giant Branch (TP-AGB) stage (e.g. Iben & Renzini 1983; 
Lattanzio & Forestini 1999; Herwig 2003). Hydrogen and helium is processed in distinct 
nuclear shells around the increasingly electron-degenerate CO-core. Periodically the He- 
shell becomes unstable due to a combination of the well-known thin-shell instability and 
partial degeneracy. The thin-shell instability is a common feature in nuclear shells around 
degenerate cores, and can be found in accreting white dwarfs (novae), in He-core white dwarfs 
as well as in accreting neutron stars (X-ray bursts). In AGB stars the He-shell flash causes a 
complicated series of mixing events, including a pulse-driven convection zone (PDCZ) which 
comprises the entire intershell layer. These mixing events eventually lead to a significant 
envelope abundance enrichment with material processed by either or both the He-shell and 
the H-shell. 

In the past most AGB stellar evolution calculations have been carried out for the solar 
metallicity, moderate metal deficiency, or zero-metallicity (e.g. Karakas et al. 2002; Siess 
et al. 2002). However, there are a number of reasons to study the evolution of AGB stars at 
very low non-zero metallicity (0.0 < Z < 0.0001). For example, recent observations of 14 N in 
metal poor stars and, in particular, in the Damped Lyman-a absorbers (Pettini et al. 2002; 
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Prochaska et al. 2002) have sparked a debate on whether the results are to be interpreted as 
a primary origin of 14 N in IMS or as a result of a peculiar (top-heavy) initial mass function 
at very low metallicity. 

Detailed abundance studies of extremely metal poor (EMP) carbon-rich stars are now 
reported in increasing numbers (e.g. Lucatello et al. 2003; Aoki et al. 2002b). Many of these 
carbon rich EMP stars are nitrogen-rich as well, and often - when it has been determined 
- these objects show elevated oxygen, sodium and s-process element abundances. A similar 
pattern with respect to C, N, and the s-process elements has been reported in mildly metal- 
poor post-AGB stars (Reddy et al. 2002; Van Winckel & Reyniers 2000). However, since 
these stars are only moderately metal poor, it is conceivable that this nitrogen abundance 
is of secondary origin, and a result of the first dredge-up (DUP). This possibility can be 
excluded for nitrogen in most C-rich EMP stars, where the N-abundance is much larger than 
the sum of the initially present C, N, and O abundances. In those cases where the binarity 
of a C-rich EMP star is established, the abundance of nitrogen and other elements may very 
well be the result of mass transfer from an AGB companion that now has evolved into a 
white dwarf. 

Intermediate-mass stars may contribute to the abundance anomalies in globular cluster 
stars in the self-pollution scenario (Cottrell & Da Costa 1981; Denissenkov et al. 1997). In 
particular, the observation of abundance anti-correlations in turn-off main-sequence stars 
(Gratton et al. 2001) has made it necessary to consider an external source for the abundance 
anomalies. Ventura et al. (2002) have proposed that the oxygen-sodium anti-correlation is 
due to the pollution of globular cluster members with material blown off TP-AGB stars. 
However, they have not explicitly included the Ne-Na and the Mg-Al H-burning cycle in 
their calculations of TP-AGB stars with Z > 0.0004. These reactions have been taken 
into account by Denissenkov & Herwig (2003), whose short study is partly based on models 
presented here in more detail. Finally, it should be noted that the TP-AGB phase of IMS 
is the dominant source of dust in early galaxies of redshifts up to 5 (Morgan & Edmunds 
2003). 

IMS during their TP-AGB evolution change their surface composition by two processes, 
the third dredge-up (DUP) and envelope burning between thermal pulses (which is also 
known as hot bottom burning, HBB). Low- mass stellar models of Pop I metallicity with 
a combination of convective overshooting and high spatial and time resolution can repro- 
duce formation of low luminosity carbon stars needed to explain the carbon star luminosity 
function in the Magellanic Clouds (Herwig et al. 1997; Herwig 2000; Mowlavi 1999). Other 
consistency checks concerning the efficiency of overshooting at the various boundaries involve 
the effects of s-process nucleosynthesis (Lugaro et al. 2003; Herwig et al. 2003). The picture 
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that has emerged from these studies of low-mass AGB stars that do not experience HBB 
is summarized in §2. In §3 the stellar evolution code and the model set is described. The 
results are presented in §4, and conclusions are given in §5. 

2. Treatment of convective boundaries 

Hydrodynamic time- and depth-dependent overshooting may be present at all convec- 
tive boundaries. An extended discussion of the motivation and background, as well as the 
technique of the exponential overshoot implementation used in this study can be found in 
Herwig (2000) and Herwig et al. (2003). The overshooting efficiency in the deep interior 
cannot be determined from hydrodynamic models ab initio due to the large range of scales 
involved. In order to construct AGB models with some predictive power, one is therefore 
forced to constrain the overshooting efficiency semi-empirically. At the bottom of the en- 
velope of TP-AGB stars, exponential overshooting increases the DUP depth at almost any 
non-zero efficiency. With an efficiency / « 0.15, the resultant 13 C pocket that forms during 
the final period of the third DUP is massive enough to reproduce the observed s-process 
overabundances (Goriely & Mowlavi 2000; Lugaro et al. 2003; Herwig et al. 2003). Models 
with overshooting at the bottom of the PDCZ predict even more efficient third DUP as well 
as a larger carbon and oxygen abundance in the intershell zone, compared with models that 
do not include overshooting at this convective boundary. The larger carbon and, in particu- 
lar, the larger oxygen abundance in the intershell is in good agreement with those observed 
in H-deficient central stars of planetary nebulae (of spectral type PG1159 and [WCj-CSPN 1 , 
Koesterke & Hamann 1997; Dreizler et al. 1996) in conjunction with the latest models on the 
evolutionary origin of these stars (Herwig et al. 1999; Herwig 2001b). These models predict 
that the PG1159 and [WC]-CSPN stars show on their surfaces the abundance distribution 
of the the intershell layer of the progenitor AGB star. In this context it is interesting to note 
that recent FUSE spectroscopic observations have established a significant Fe-deficiency in 
some hot H-deficient central stars (Miksa et al. 2002). Further observations may substantiate 
the hypothesis that this Fe-depletion is a result of the s-process in the intershell during the 
progenitor AGB evolution. This would support the claim that hot H-deficient central stars 
of PN display directly the intershell material of their AGB progenitors. 

However, in the s-process models, the larger carbon intershell abundance in models with 
overshooting leads to a larger neutron exposure in the s-process layer during the interpulse 
phase (Lugaro et al. 2003). In these models the s-process abundance distribution is too top- 



Central stars of planetary nebulae with Wolf-Rayet type spectrum. 
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heavy, or in other words the [hs/ls] index is larger than observed. This could be compensated 
by some weak mixing of the s-process layer during the interpulse phase (e.g. as a result of 
slow rotation, Herwig et al. 2003). 

Models with exponential overshooting at the bottom of the PDCZ predict larger tem- 
peratures in the second s-process production location, fuelled by the 22 Ne neutron source. 
Some branchings like the 96 Zr/ 94 Zr ratio are sensitively dependent on the temperature at 
the bottom of the PDCZ. At the same time, such isotopic ratios are experimentally known 
to high precision from mainstream pre-solar SiC grains, steming from low-mass AGB stars 
(Zinner 1998; Nicolussi et al. 1998). This puts perhaps the most stringent constraint on the 
maximum possible overshooting efficiency Although a more detailed study in this respect is 
still pending, preliminary calculations by Lugaro et al. (2003) indicate that low mass models 
with /pdcz = 0.008 may be in agreement with the temperature-dependent branchings, and 
at the same time retain the enhanced oxygen intershell abundance indicated by the PG1159 
and [WC]-CSPN stars. 



3. Stellar evolution code and model calculations 

The models in this study were computed with the full ID stellar evolution code of Herwig 
(2000), including the exponential overshooting model described there. An earlier version of 
this code has been used by Blocker (1995) to study the evolution of Pop. I AGB stars as 
a function of initial mass. As in most other AGB models the Schwarzschild critierion for 
convection is used. The stabilizing effect of a /z-gradient as reflected by the Ledoux criterion is 
thus neglected. Using the Ledoux criterion in AGB stars would make it even more difficult 
to obtain the third dredge-up. In addition, the exponential overshooting included here, 
would in this situation likely dominate mixing through semiconvection. In order to meet 
the high numerical resolution requirements of third DUP models, the adaptive time step 
and grid size routine was improved. It can now deal with the hot DUP burning episodes 
described below. In particular, it is ensured that the time-step and the grid resolution 
correspond to each other in the region near the bottom of the convective envelope during the 
DUP phase. In this situation, the Lagrangian advance rate of the bottom of the convective 
envelope mcE implies a relation between the grid resolution and the time step. Numerical 
experiments indicate an appropriate Lagrangian scale in the third DUP overshooting zone 
to be / m = 10~ 6 M Q , which should be well sampled in space and time. This means that the 
typical grid spacing in the DUP zone is < 10~ 7 M Q and time steps should not be larger than 
^m/^CE- In practice, this means that a single third DUP episode of several 10~ 3 M requires 
computing several thousand models. In the case where nuclear burning takes place during 
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the DUP phase - as described below - the time step is further limited by the fluctuations of 
the hydrogen-burning luminosity L H . 

The opacities are from Iglesias & Rogers (1996, OPAL) supplemented with low temper- 
ature opacities by Alexander & Ferguson (1994). The mixing-length parameter has been set 
to «mlt = 1-7, as calibrated by fitting a solar model. The nuclear network equations and the 
time-dependent mixing equations are solved simultaneously, using a fully implicit, and itera- 
tive scheme. For this study an automatic sub-time step algorithm has been implemented into 
the Newton-Raphson scheme. If the convergence criteria are not meet within five to eight 
iterations the time step is replaced by 10 times smaller sub-time step, followed by a second 
sub-time step that takes care of the remaining 90% of the original time step. This procedure 
can be applied repeatedly, and guarantees that a solution with the pre-defined precision can 
be found even in fairly violent burning and mixing events, as encountered in massive AGB 
models of very low metallicity. All relevant charged particle reactions are included up to and 
including the Mg-Al cycle, as well as neutron captures. Reaction rates are mainly based on 
Caughlan & Fowler (1988). Updates from El Eid & Champagne (1995) and from Angulo 
et al. (1999) for the important CN-cycle p-capture rates have been considered. 

The benchmark model sequence (E79-D4, see Table 1) has an initial mass of 5M and 
a metallicity of Z — 0.0001 ([Fe/H] m \og(Z/Z Q ) = —2.3). It is evolved from the pre-main 
sequence through all evolutionary phases up to the AGB. For the initial helium mass-fraction 
abundance we adopt Y = 0.23025 according to 

Ay 

Y = Y Q + —x(Z-Z G ), 

where Y & = 0.28 is the value adopted for the solar initial helium abundance (Grevesse & Noels 
1993; Grevesse & Sauval 1998), and Z Q = 0.02 has been adopted for the solar metallicity. 
AY/AZ is the ratio of fresh helium to metals supplied to the interstellar medium by stars. 
The value of 2.5 used in this study is consistent with many observational constraints (e.g. 
Pagel & Portinari 1998, and references therein). 

Many stars of low metallicity may, in fact, form from gas clouds that have a non-solar 
abundance distribution. For example, Carney (1996) concluded in his review of observational 
data that globular cluster stars with [Fe/H] ~ —0.6 are overabundant in the a-elements by a 
factor of about 2. However, some halo stars have been reported with solar-scaled a-element 
abundances (Nissen & Schuster 1997). For this study we have adopted a scaled solar metal 
distribution according to Grevesse & Noels (1993) and Anders & Grevesse (1989). Most of 
the results discussed here are related to the production and processing of primary C and O, 
which dominate the surface abudance evolution once DUP starts. 

Mass loss has not been considered for the 5 M cases, on which most of this investigation 
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is based. This choice was motivated by the fact that model tracks with mass loss show less 
DUP (Karakas et al. 2002). Besides, the amount of mass loss that occurs in very low 
metallicity stars is uncertain. Due to the lack of mass loss the calculations do not have a 
natural end due to the loss of the stellar envelope. At least one sequence is followed through 
15 TPs which is sufficient to reach the asymptotic regime in which changes of many quantities 
are smaller and more regular than during the initial TP-AGB phase. 

Up to the first TP on the AGB, exponential overshooting with / = 0.016 is applied 
at all convective boundaries. This is the same value as used by Herwig (2000), and it 
was originally obtained by calibrating core-overshooting to reproduce the observed width of 
the main-sequence. For the TP-AGB benchmark case (sequence E79-D4) a non-zero but 
negligible overshooting efficiency 2 of /pdcz = 0.002 at the bottom of the PDCZ and the 
previously used value of /ce = 0.016 at the bottom of the convective envelope are applied 
at all times - during the HBB interpulse phase as well as during the DUP phase. In our 
opinion, this choice represents a minimum overshooting situation. AGB models with such 
overshooting parameters would not generate enough 13 C for the s-process production during 
the intershell, and they would not be able to reproduce the observed abundances of PG1159 
stars and [WC]-CSPN (§2). For this reason, the possibility of larger overshooting efficiency 
is explored by a number of test computations (Table 1, §4.5). Rotation and the possible 
effects of magnetic fields are not included. 



4. Results 

4.1. Pre- AGB evolution 

The 5 M Q pre-main sequence model is evolved through all phases up to the giant stage. 
The evolution track is in good agreement with calculations by Girardi et al. (1996) (Fig. 1). 
The small remaining differences can be attributed to differences in overshooting efficiency, 
mixing-length parameter, and other details of the physical inputs and numerical techniques. 
The core H- and He-burning times - in particular the sum of both - are in reasonable 
agreement as well. This model has th = 8.64 x 10 7 yr and rn e = 1-21 x 10 7 compared to 
r H = 9.00 x 10 7 yr and r Hc = 9.13 x 10 6 yr for the model by Girardi et al. (1996). Initially 
the mass of the main-sequence core convection zone is 1.33 M , and decreases to ~ 0.45 M 
shortly before H-burning in the core ceases. The maximum extent of the convective He- 



2 This small amount of overshooting has very little effect on the He- flash peak luminosities and intershell 
abundances, but behaves numerically more benign with the simultaneous solution of time-dependent mixing 
and nuclear burning. 
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burning core is 0.55 M and one major breathing pulse has been found. At the end of 
the core He-burning phase, the H-shell is located at m r = 1.47 M . The second DUP starts 
4.25 x 10 6 yr after the end of the He-core burning when the bottom of the convective envelope 
reaches into the H-free core. This mixing episode ends after 63000 yr when the bottom of 
the envelope has reached m r = 0.9615 M Q . 

Intermediate-mass stars at this low metallicity do not have a first giant branch phase 
(Girardi et al. 1996) and thus the first envelope abundance alteration occurs as a result of 
the second DUP. Material processed mainly by H-shell burning is brought to the surface. 
4 He is enriched by about 50% to X( 4 He) = 0.335 and hydrogen is reduced accordingly to 
X(H) = 0.65. Compared to the very efficient production of 14 N by the combined action of 
HBB and the third DUP later during the TP-AGB evolution, the second DUP of 14 N with 
an enhancement by ~ 50% can be considered unimportant. 23 Na is enhanced by 0.75dex 
from the conversion of the initial abundance of 22 Ne and some 20 Ne (see also §4.4.3). 16 is 
depleted by O.ldex due to partial ON cycling. Magnesium isotopes are changed by less than 
0.05, 0.1 and O.Oldex for mass numbers 24, 25 and 26 in the second DUP. The first He-shell 
flash occurs 12200 yr after the end of the second DUP. At this time the H-free core has a 
mass of 0.9630 M m . 



4.2. Envelope burning and third dredge-up 

The entire chemical enrichment is strongly dependent on the occurrence of the third 
DUP. The 5 M Q TP-AGB model sequence shows efficient third DUP after the fourth thermal 
pulse (Fig. 2). The DUP mass is almost AMqup = 3 x 10~ 3 M Q , comparable to the DUP in 
a low-mass AGB star with A = 0.5. This DUP parameter is defined as A = AMh/AMqup, 
where AMh is the core mass growth due to H-burning during the interpulse phase. All 
following thermal pulses have a very efficient third DUP as well. The fifth thermal pulse is 
an exception because the time-step control algorithm failed. As a result, the DUP efficiency 
is too small (Fig. 2). For all subsequent TPs the DUP parameter always exceeds unity, 
and the H-free core M H effectively shrinks slightly. This is also true for the He-free core 
Mn e . The evolution of the He- and H-burning luminosity, as well as the stellar radius and 
the stellar luminosity, of the 5 M model is shown in Fig. 3 and 4. Note how the surface 
properties as well as the evolution of the shells, change at the onset of deep dredge-up. The 
evolution of shell luminosities and surface parameters differs strongly from the situation in 
more metal-rich and/or less massive models (Lattanzio 1992; Blocker 1995; Herwig 2000). 

A more detailed presentation of the sixth thermal pulse is given in Fig. 5. The evolution 
of the bottom of the envelope convection zone is very smooth. In various test sequences the 
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time evolution of this convective boundary showed discontinuous step-like advances into the 
core. In these tests the numerical resolution was insufficient. Such features cannot be found 
during the DUP after the sixth or any subsequent thermal pulse of the benchmark 5 M Q 
sequence. In this model sequence the third DUP, in particular, with respect to the Lagrange 
location of the He-burning shell (Mn e ) and the PDCZ, is very deep. The envelope convection 
not only engulfs the entire interpulse zone, but it even penetrates into the He-free core. This 
third DUP behavior is not present in more metal-rich massive AGB models or in low-mass 
AGB models. A decreasing M He has previously been reported by Herwig (2000, Table 1) for 
the TP- AGB evolution of a M ini = 4 M sequence with Z = 0.02. However, for that sequence, 
efficient overshooting (/pdcz = 0.016) at the bottom of the PDCZ has been assumed, which 
leads to a deep penetration of the PDCZ into the C/O core beneath the He-shell. This is 
not the case in the present model, which has been computed with negligible overshooting at 
the He-shell flash boundary (/pdcz — 0.002). The bottom of the PDCZ does not penetrate 
significantly into the He-shell, and consequently the intershell abundances are not enhanced 
with carbon and oxygen, as has been found in models with efficient PDCZ overshooting. 
The mass fractions in the PDCZ of the sixth TP are ( 4 He/ 12 C/ 16 0) = (0.854/0.143/0.002) 
(see below for the tenth TP). 

In all successive thermal pulses deep DUP reaching below the He-shell is present, al- 
though the mass dredged- up from below the He-shell decreases somewhat (Fig. 2). The 
events during the tenth thermal pulse are shown in detail in Fig. 6. The four abundance 
profile panels (i-iv) show the situation at the successive times indicated in the top panel. 
Panel (i) shows the situation shortly before the onset of the He-shell flash. The dip in the 
16 O profile at location (A) is a signature of the H-shell. In the convective envelope (light 
shaded region) 14 N is the most abundant species of the CNO elements, indicating efficient 
HBB. At location (B) the large abundance of 14 N reveals the nature of this material as H- 
shell ashes. Toward location (C) 12 C, 16 0, and as well 22 Ne gradually increase in abundance, 
while 14 N decreases as a result of a-capture reactions. Location (C) harbors He-shell ashes 
which constitute the top layer of the C/O core. 

Shortly after t = 49.51 x 10 3 yr, the He-shell flash drives the PDCZ. Panel (ii) shows the 
abundances at the peak extension of the PDCZ (darker shading). Just above m r = 0.9615 M 
the tiny He-buffer and the signatures of the H-shell can be seen. Thus, the PDCZ does not 
reach into the H-rich envelope and no protons are ingested into the convection zone. Such an 
occurrence has been suggested to be an alternative site for the s-process in extremely metal 
poor or zero-metallicity TP- AGB stars (Aoki et al. 2001). If such a mixing event occurs in 
intermediate mass stars it should be restricted to metallicities lower than that studied here 
(Z=0.0001). 
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During the pulse 14 N is synthesized into 22 Ne, and to some extent into 25 Mg. The latter 
is the result of the 22 Ne(a, n) 25 Mg reaction, which releases neutrons for the s-process. In 
addition, 12 C and 16 are now present in the intershell with a larger abundance than before 
the He-flash (( 4 He/ 12 C/ 16 0) = (0.782/0.198/0.017)). The 16 abundance now exceeds 1% by 
mass, a value which is very similar to predictions of more metal-rich models obtained without 
overshooting (Schonberner 1979). We discuss in §4.5 how the abundance of this primary 
oxygen depends on the overshooting. However, the presence of 1 — 2% of oxygen by mass in 
the intershell of AGB stars after a few initial TPs should be considered as a minimum at any 
metallicity. The third DUP will mix this primary oxygen into the envelope. In stars of solar 
or only mildly sub-solar metallicity, the presence of primary oxygen might be unnoticeable 
observationally, partly because of the inherent difficulties in determining oxygen abundances, 
and partly because the initial oxygen abundance is very close or even equal to the oxygen 
abundance in the intershell. However, this obviously changes in TP-AGB stars of very low 
metallicity. If the initial envelope abundance has extremely low CNO abundance, then the 
DUP of material containing X( 16 0) = 0.02 should be noticeable. In fact, the zero-metallicity 
models of both Chieffi et al. (2001) and Siess et al. (2002) report a significant fraction of 
oxygen in the dredged-up material that is especially obvious in the less massive cases without 
very efficient HBB. From the intershell abundances reported here, one can immediately make 
some simple predictions. For example, if a hypothetical extremely metal poor star has a C- 
overabundance of [C/Fe]=2.5, and if this C-abundance comes from the heavy pollution by 
a TP-AGB star that has not experienced very efficient HBB (which would have destroyed 
some O) then the prediction for oxygen would be [O/Fe] > 1.5. This is so simply because 
the predicted X c /X ratio in the dredged-up material is < 10. 

In panel (iii) the situation during the third DUP following the He-flash is shown. In 
location (D) the bottom of the envelope proceeds inward, bringing to the surface 12 C, 16 0, 
22 Ne, and other species. This 22 Ne is the seed for 23 Na production by HBB. In location 
(E), panel (iii) a small 25 Mg pocket has formed from the a-capture reaction on 22 Ne, that 
releases neutrons for the s-process. The neutrons released during the production of this 
25 Mg pocket are in addition to the neutrons from the 22 Ne neutron source considered in the 
present s-process models. These neutrons are released during a high-T radiative 22 Ne burning 
phase (along the dashed line in the top panel of Fig. 6) between t = 49.515 x 10~ 3 yr and 
t = 49.53 x 10~ 3 yr. The temperature during this 15 yr period of time decreases from the peak 
PDCZ value of logT = 8.55 to logT = 8.3, and for 10 out of the 15 yr, logT > 8.4. In more 
metal rich or low mass TP-AGB models, the DUP never reaches this layer. Consequently, 
present s-process models do not need to consider these neutrons, which usually are buried 
in the core. However, in this massive TP-AGB model of very low metallicity, the DUP picks 
this material up as it reaches below the He-shell. This is shown in panel (iv), where the 
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situation at the end of the third DUP is shown. 

Two interesting things can be seen in panel (iv). Firstly, because the DUP now reaches 
below the He-shell, additional primary oxygen is brought to the envelope in significant 
amounts. Thus, oxygen is dredged-up efficiently, despite the fact that practically very 
inefficient overshooting has been assumed at the bottom of the PDCZ. Secondly, the 14 N 
abundance in the bottom region of the convective envelope exceeds the C and O abundance. 
This means that H-burning already takes place at the bottom of the envelope. HBB has 
already started and processes C and O dredged-up from the core in situ. 

This dredge-up is hot as opposed to the situation in less massive or more metal rich TP- 
AGB stars. In sequence E79-D4 the peak H-luminosity during the efficient hot DUP events 
is in the range 1.4 — 1.8 x 10 6 L and is not changing in a systematic way from pulse to pulse. 
The large H-luminosity leads to a peak in the surface luminosity (Fig. 4). At low metallicity 
the envelope reaches faster into deeper and hotter layers than at higher metallicity. The 
temperature in the overshooting layer during the DUP remains below the value obtained 
in the H-shell during subsequent HBB. However, the 12 C abundance is much higher in the 
intershell then in the H-shell. Therefore H-burning is not limited by the slow p-capture of 
14 N, but only by the faster 12 C(p,7) 13 N reaction. As a result nuclear energy production by 
H-burning in the overshooting layer during DUP is several orders of magnitude larger then 
in the H-shell during HBB, and in particular vigorous in the exponential overshooting zone. 
However, even within the convective envelope the H-burning luminosity is still significant. 
In this innermost region of the convective envelope, the abundance of 12 C and 14 N are deter- 
mined by both mixing as well as p-capture nucleosynthesis. Due to the coupled treatment of 
time-dependent mixing and nucleosynthesis, the models predict a continuous decrease from 
the large intershell values to the small envelope abundance. In the transition zone at the very 
bottom of the envelope convection, the H-burning luminosity decreases gradually outward 
from the peak value in the overshooting layer. For large core masses like for the benchmark 
5 M Q case E79 this additional luminosity increases the DUP efficiency. Since the radiative 
temperature gradient is proportional to the luminosity (V ra d ~ /) efficient H-burning during 
the DUP increases the convective instability, and models including this effect show more effi- 
cient envelope enrichment at each TP. In § 4.5 the dependence of this effect on the overshoot 
parameter is discussed. 

The H-burning during the DUP has been noted before for Z=0 TP-AGB models with 
a small amount of overshooting by Chieffi et al. (2001). Their H-burning luminosities are 
even similar. However, only one or two pulses with this hot DUP have been followed due to 
the numerically demanding nature of this phase. 

In addition to the Z = 0.0001 models, a comparison sequence with Z = 10~ 5 has 
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been computed. Overall the properties are very similar. In particular, the deep DUP, the 
hydrogen burning during the DUP phase, and the efficient 14 N production through HBB is 
found in the intermediate-mass stars of the lowest metallicities. In the extremely metal poor 
cases, the envelope CNO abundance is dominated by the primary component from the third 
DUP. It is this envelope CNO abundance which determines to a large extent the evolution 
of envelope burning AGB stars, including their yields. Since DUP is so efficient and sets in 
very early, stars with even lower metallicity than Z = 0.0001 still behave in a very similar 
way. 

4.3. s-process considerations 

The main nuclear production site of the s-process elements is usually associated with 
low-mass TP-AGB stars (Gallino et al. 1998; Busso et al. 1999; Goriely & Mowlavi 2000; 
Busso et al. 2001; Lugaro et al. 2003; Herwig et al. 2003). This, however, does not imply 
that massive AGB stars do not synthesize the s-process elements at all. Very little is known 
theoretically about the s-process at very low metallicity. IMS of this metallicity may have 
a distinctive s-process signature that can be used as a diagnostic tool for the interpretation 
of observations. Two reactions can release significant amounts of neutrons in conditions 
encountered in TP-AGB stars. The 13 C(o;,n) 16 reaction is efficient for logT > 7.95 while 
the 22 Ne(o;, n) 25 Mg reactions requires logT > 8.45 to release significant amounts of neutrons. 

According to current models of the s-process in low-mass stars 13 C releases neutrons 
under radiative conditions during the interpulse phase and provides most of the neutrons 
in the majority of s-process enriched stars. It is also believed that this neutron source is 
responsible for s-process isotopic signatures found in pre-solar meteoritic SiC grains. For this 
mechanism, protons have to be partially mixed with 12 C at the intershell-envelope interface 
at or immediately after the end of the third DUP. This partial mixing zone will heat up as 
the intershell layer contracts to the pre-flash configuration. Then 13 C forms which at a later 
phase during the interpulse will capture a-particles and release neutrons. This scenario has 
been very successfully applied in conjunction with stellar evolution models that do assume 
an ad-hoc H-profile for the partial mixing zone. As discussed in § 2 recent studies of the 
effects of rotation and/or overshooting as a physical mechanism for the required mixing have 
revealed a more complicated picture that is not well understood yet. 

In the very low metallicity IMS with deep DUP, a partial mixing zone forms at the end 
of the third DUP as a result of the overshooting algorithm. The first difference compared 
to the low-mass models is that the 13 C pocket forms immediately. As shown already in the 
previous section, the temperature at the bottom of the convective envelope is high enough 
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for H-burning during the actual DUP period. The narrow 13 C pocket at the end of the third 
DUP episode of the tenth TP (corresponding to panel (5) in Fig. 6) is shown in Fig. 7. The 
14 N pocket which acts as a neutron poison is also evident. Since the DUP has reached below 
the He-shell, the mole fraction of 4 He (for the subsequent activation of the 13 C(o;,n) 16 
reaction) is only a fraction of the 13 C abundance in the pocket. For this reason, the 13 C 
neutron source may be not or only inefficiently available in these massive TP-AGB models. 

The (a, n) reaction of 22 Ne releases neutrons in the He-shell during the He-flash. This 
neutron source is activated very efficiently. 14 N from which 22 Ne is generated is the result 
of CNO cycling on previously dredged- up primary C and O. The neutron source efficiency 
therefore depends on the peak temperature that is reached as well as the cumulative efficiency 
of previous DUP episodes. In this model the third DUP is very efficient and therefore the 
14 N abundance in the ashes of the He-shell is rather large. This enhances the neutron release 
in the PDCZ, and this effect is obviously more pronounced at lower metallicity. A second 
contribution comes from radiative 22 Ne burning already mentioned in §4.2. Following the 
peak of the He-shell flash the temperature in the He-shell gradually decreases. At this time 
the actual He-shell is no longer convectively unstable. Neutron captures may still be possible, 
but now under radiative conditions. Later this layer will be picked up by the third DUP. This 
second contribution will likely not exceed 10% of the s-process synthesis in the PDCZ, but 
may be important for details at certain branchings. Finally, it may be worth considering the 
effect of 13 C left behind in the H-shell ashes. Since the primary CNO abundance can be quite 
large due to efficient DUP, the ratio of 13 C to 56 Fe seed will be much more favorable than 
in more metal-rich stars. On the other hand, due to the large 14 N/ 13 C ratio in the H-shell 
ashes, the efficiency of the trans-iron element production will be limited. A quantitative 
analysis of the s-process in these models will be presented elsewhere. 



4.4. Surface abundance evolution 

4-4-1- H and He 

During the TP-AGB phase the envelope abundance is altered by the combined effect 
of DUP and envelope burning. After 14 thermal pulses the envelope hydrogen abundance 
of the E79-D4 sequence has decreased from initially A(H) = 0.665 to 0.654. Initially the 
H-abundance changes mainly because of the third DUP. However, after the eighth thermal 
pulse, the H-destruction due to HBB dominates and during the following TPs the H-depletion 
follows a roughly linear decrease with time that can be approximated by 

X(H) = -1.3 x 10~ 4 x t 3 + 0.6669 
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with i 3 = t/1000. Correspondingly the 4 He abundance increases from X( 4 He) = 0.335 to 
0.343 and follows approximately the relation 

X( 4 He) = 1.15 x 10" 4 x t 3 + 0.3327. 

4-4-2- C, and primary N 

The combined CNO abundance increases steadily and follows the relation 

X(C + N + O) = 2.1 x 10~ 5 x t 3 + 4.0 x 10 4 . 

Thus, at each thermal pulse, the mass fraction of CNO material in the envelope increases 
by AX R3 2 x 10~ 4 . The relative abundance of the CNO isotopes changes periodically with 
the pulse phase (Fig. 9). After the onset of efficient third DUP, the temperature during the 
interpulse phase at the bottom of the convective envelope (HBB temperature) quickly settles 
in the range 9.1 — 9.6 x 10 7 K. Primary C and O dredged-up from the intershell is cycled 
into 14 N. Part of that cycling occurs immediately during the DUP, which can be seen from 
the increase of 13 C together with 12 C and 16 at the DUP time (§4.2). During this short 
H-burning DUP phase, the 12 C/ 13 C ratio rises to values in the range 6 — 7, which is much 
larger than the CNO cycle equilibrium value of ~ 3.5 that is obtained in the envelope in the 
second half of the interpulse cycle due to HBB. 

Even though the HBB temperature is high enough to burn oxygen, the effect of oxygen 
DUP is clearly noticeable. In this model, third DUP of oxygen prevents an efficient 16 O 
depletion by HBB. The oxygen isotopic ratios evolve towards 16 0/ 17 > 100 during the 
interpulse phase, but larger values (120 — 180) are present in the envelope during the earlier 
TPs, and in particular during the H-burning DUP events. The 16 0/ 18 > 2 x 10 6 at all 
times. 

The most important aspect of the CNO abundance evolution is the significant primary 
14 N production. After the 14 TPs computed here without mass loss, the 14 N abundance is 
twice the solar 14 N abundance. The envelope, that will be lost eventually, of this model star 
contains at this point approximately 8 • 10~ 3 M Q of primary 14 N. 

4-4-3. Oxygen and Sodium 

The elemental abundance evolution of sodium and oxygen is correlated or stationary 
in all models reported here. Oxygen in the envelope increases because of dredged-up, and 
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decreases because of HBB. In this model these two effects balance each other. Sodium forms 
early on in the pulse cycle, and is then destroyed again. This concordance can easily be 
explained. 22 Ne - the seed for 23 Na - follows indirectly the dredge-up of primary C and 
O from the intershell, as it is made from the 14 N that results from H-shell burning of the 
dredged-up C and O material. If the envelope burning temperature is large enough for the 
22 Ne(p, 7) 23 Na reaction, then the production of 16 and 23 Na are correlated. The destruction, 
on the other hand, is also correlated because the p-capture rates of both isotopes are very 
similar (Angulo et al. 1999). 

Denissenkov & Herwig (2003) have discussed the implications of these processes in 
massive AGB stars for the O-Na anti-correlation in globular clusters. They have argued 
that low-metallicity IMS are probably not directly responsible for the observed O-Na anti- 
correlation. An important detail in their argument is the recurrent DUP of oxygen (§4.2), 
which would require larger HBB temperatures for effective oxygen depletion. At such large 
temperatures any sodium produced by p-capture on 22 Ne would be destroyed as well. In a star 
with lower HBB temperatures, simultaneous oxygen and sodium enhancement is predicted. 
In any case, AGB stellar models predict an O-Na correlation, or no correlation, rather than 
an anti-correlation. In Fig. 10 the sodium and oxygen abundances the globular cluster stars 
and the model prediction from this work are compared. The globular cluster stars show a 
varying degree of oxygen depletion that is associated with a sodium enhancement. 

The offset in sodium abundance between the model abundances and the globular clus- 
ter stars can be fully attributed to the choice of initial composition for the models. The 
TP-AGB initial sodium abundance of the models, which correspond to [Na/Fe] ~ 0.6 and 
[O/Fe] ~ —0.2 (Fig. 10), are the result of the 2 nd DUP, during which the envelope convection 
reaches into the ashes left behind by the H-shell. In the H-shell, primordial 22 Ne is partially 
transformed into 23 Na. The initial sodium abundance of the TP-AGB model corresponds 
to an initial scaled-solar 22 Ne abundance. While no observational information on the noble 
gas neon is available for low metallicities, theoretical models of galactic chemical evolution 
agree that [Ne/Fe] is larger at lower [Fe/H]. At [Fe/H]= —2.5 the predictions range from 
[Ne/Fe]= +0.2 (Timmes et al. 1995) to [Ne/Fe] = +0.7 (Alibes et al. 2001). However, this is 
predominantly 20 Ne, since 22 Ne is a secondary isotope and follows the initial abundance of 
14 N in massive stars (Woosley & Weaver 1995). Unless the parent cloud of the AGB stars 
modeled here has been polluted by an even earlier generation of massive AGB stars, its 22 Ne 
abundance has probably been negligible. In that case, the initial sodium abundance of the 
Z=0.0001 massive TP-AGB stellar models is the same as the initial sodium abundance at the 
time of star formation. In order to explain the O-Na anti-corelation in globular clusters by 
simultaneous Na production and O destruction it is often assumed that [Na/Fe] < initially 
at low metallicity. This is in agreement with theoretical models of the galactic chemical 
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evolution (Timmes et al. 1995; Alibes et al. 2001), although a systematically lower [Na/Fe] 
at lower [Fe/H] is not obvious from the analysis of field giants by Gratton et al. (2000). 
However, they report a scatter of ~ ±0.3 around [Na/Fe] = for —2.5 < [Fe/H] < —0.5. 

4-4-4- Magnesium 

The magnesium isotopic ratios can be modified by three processes in AGB stars: (a) 
a-captures on 22 Ne in the He-shell (flash), which can produce both 25 Mg and 26 Mg, (b) 
p-captures on 24 Mg and 25 Mg during the envelope burning phase, and (c) by n-captures on 
24 Mg and 25 Mg as well as on 26 Al. In Denissenkov & Herwig (2003) the TP-AGB evolution 
dominated by high HBB temperature is discussed, during which the initially available 24 Mg 
burns mainly into 25 Mg. In cases of lower HBB temperature (either because of the inclusion 
of mass loss or because of a smaller mass), and if a complete treatment of neutron capture 
reactions is taken into account, the other two effects operating in the PDCZ become more 
important for the model predictions. 

4.5. Dependence on the overshooting parameter 

Three additional tracks with different assumptions on the overshooting efficiency have 
been calculated (Table 1). Comparison of sequence D4 and D9 reveals the influence of the 
overshooting parameter at the bottom of the PDCZ. For sequence D9 a larger overshooting 
efficiency (/pdcz = 0.016) than for sequence D4 has been assumed. These tests confirm 
the trends found for the 3M Q , Z=0.02 stellar model with intershell overshooting by Herwig 
(2000). In comparison to the D4 sequence, the D9 sequence shows a 

1. larger Lne-flash peak luminosity by roughly a factor 1.5 — 2 

2. longer interpulse period by ~ 20% 

3. slightly smaller DUP H-burning luminosity, and similar quiescent interpulse H-shell 
luminosity 

4. larger third DUP, e.g. at the sixth TP of the D9 sequence at t — 1 — 35384 yr (t as in 
Fig. 2) the DUP parameter is A = 1.46 compared to A = 1.18 at the eighth TP of the 
D4 sequence at t — to = 32644 yr; cumulatively this deeper DUP leads to a somewhat 
stronger core mass reduction than in the D4 case; the core mass at the end of the sixth 
TP of the D9 sequence is m T = 0.95539 M , and 



-17- 



5. a different He/C/O intershell abundance composition; C and O are enhanced at the 
expense of He; in the PDCZ of TP6 the mass fractions are (He/C/O) = (0.62/0.26/0.18). 

Despite these differences, overshooting at the bottom of the PDCZ does not change the 
envelope abundance evolution qualitatively. In both cases, with and without overshooting 
at the bottom of the PDCZ, the DUP is efficient, and substantial amounts of primary C and 
O will be dredged-up. Sequences of intermediate-mass and massive AGB models generally 
tend to converge to A = 1 whether they start initially with A > 1 or A < 1 (Karakas et al. 
2002). This phenomenon is already noticeable in Fig. 2. Therefore, the difference in DUP 
efficiency is expected to decrease in later pulses. The different 16 and 12 C abundances in 
the PDCZ lead to some differences in the surface evolution of the CNO elements. For a 
comparison, the abundance evolution of the D4 and D9 case are shown in Fig. 11 in the 
logaritmic square bracket units used to present abundance observations. Deeper DUP and 
a larger 12 C-abundance in the PDCZ cause a larger enrichment in primary 14 N in the D9 
case compared to the D4 case. A difference of the 16 abundance of a few tenths of a 
dex is present. If many more TPs occur in this type of stars, the difference of the oxygen 
abundance in the two cases would likely increase. However, overall the uncertainty in the 
model predictions of massive AGB stars do not appear to be strongly dependent on the 
overshooting efficiency at the bottom of the PDCZ. 

Next, the influence of the exponential overshooting efficiency at the bottom of the con- 
vective envelope during the DUP phase is analyzed. A number of tests (among them E79-D10 
and Dll with f C E > 0.032) have been followed for several ten thousand models with small 
time steps during the hot DUP phase. In some cases the mixing speed has been altered 
too. The dredge-up rate was almost constant in all cases (dm^ ce /dt ~ 5 • 1O~ 5 M / yr for 
E79-D10, with m^ ce is the Lagrange coordinate of the bottom of the convective envelope), 
and depends somewhat on both the assumed speed and depth of the exponential overshoot- 
ing. The surface abundance evolution during the ongoing hot DUP after the fifth TP (e.g. 
sequence E79-D10 which has been followed longest, Fig. 12) is similar to the long-term evo- 
lution over many TPs with moderate DUP (Fig. 11). However, in the latter case the N/O 
and the N/C ratios are larger and the C/O ratio is smaller. 

When the sequence was stopped the dredged-up mass was AMdup = 1-6 • 1O~ 2 M , 
which is about eight times the mass contained in the He-shell flash convection zone. During 
the hot DUP the stellar luminosity is about logL/L Q = 5, again somewhat dependent on 
the details of the overshooting. With the large envelope CNO abundance enhancement the 
mass loss may not be less efficient than in more metal-rich AGB stars. Assuming the mass 
loss rate of Blocker (1995) with an efficiency parameter 7/bl = 0.1 the mass loss during this 
phase is of the order M — 2 ■ 1O _3 M / yr. Together with an envelope mass of about 4M 
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this implies an upper limit for this evolution phase of toup < 2000 yr, and with the dredge-up 
rate the maximum mass of dredged- up material is AM DUP ~ 0.1 M . 

It is emphasized that these findings are the result of test calculations intended to explore 
the uncertainties of the model predictions. Based on the current models, and without good 
constraints on overshooting in the presence of nuclear burning, we can not exclude that IMS 
of very low metallicity end their lives in a high mass loss hot-DUP phase, possibly after a 
rather early thermal pulse. 

5. Conclusion 

The evolution of massive AGB stars (mainly of initial mass 5 M and very low metallic- 
ity, Z = 0.0001) has been studied in detail. Nucleosynthesis and convective mixing, including 
the effect of exponential overshooting has been considered. In contrast to less massive cases 
and more metal rich cases this model predicts hot DUP, the efficient nuclear burning of 
hydrogen during the third DUP phase which can greatly increase the DUP efficiency The 
additionally released energy feeds back into the structure evolution, leading to deeper DUP 
than observed in models that do not show this effect (either because strictly no mixing 
beyond the Schwarzschild boundary of convection is allowed or because the mass and metal- 
licity are not in the appropriate range). In the benchmark case a DUP penetration through 
and beneath the He-shell is found even with a rather small overshooting efficiency. For larger 
masses or larger overshoot efficiencies, the DUP is deeper and may eventually terminate the 
AGB evolution. 

Our model predictions of O and Na in massive low-Z TP-AGB models is not in agree- 
ment with the trend seen in the globular clusters, and this casts doubt on the simple self- 
pollution scenario of globular clusters. On the other hand AGB stars have been frequently 
associated with the significant overabundance of C and N in some EMP stars. EMP bi- 
naries CS 29497-030 (Sivarani et al. 2003) and LP 625-44 (Aoki et al. 2002a) show oxygen 
and sodium enhancements similar in magnitude to the model predictions presented here. 
However, the majority of C-rich EMP (including these two) stars have a larger C than N 
enhancement. Such an abundance pattern is incompatible with HBB in massive AGB stars, 
in which at all metallicities C and O is cycled efficiently into N (Fig. 11 and 12). 

The EMP star HE 0024-2523 (Lucatello et al. 2003) is a post-common envelope binary 
with a period of only a few days. Our models suggest that the most likely time to initiate 
the common envelope evolution should be at the time of the first thermal pulse with efficient 
DUP. In the 5M sequence this is the fifth TP. During this phase the star expands tem- 
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porarily very significantly (Fig. 4). The radius peaks occur at all subsequent thermal pulses, 
but the maximum radius during these peaks is not significantly exceeding the initial radius 
peak. Thus, if the common envelope evolution has not been initiated during the preceeding 
RGB or early AGB evolution, the most likely time is at the first TP with efficient DUP. We 
do not model the common envelope evolution in this paper. However, it is interesting to 
note that the C/O ratio in the intershell is about 100 during the first TP with deep DUP 
§4.2. During the following TPs the intershell C/O ratio decreases to ~ 15 at the tenth TP. 
With this C/O ratio in the intershell abundance the observed C/O ratio in HE 0024-2523 of 
100 can not be reproduced no matter how peculiar the burning and mixing induced by the 
common-envelope evolution may be. However, HE 0024-2523 is highly s-process enhanced 
and it is not clear how this can be explained without a succession of neutron exposures in 
many subsequent TP cycles. This example shows that the s-process in AGB stars of very 
low metallicity is not very well understood yet. 

Our calculations confirm the previously found trend that both a lower metallicity as 
well as a larger core mass increase the DUP efficiency in terms of the DUP parameter A 
(e.g. Boothroyd & Sackmann 1988). Marigo (2001) and Marigo et al. (1999) found semi- 
empirically that low-mass and moderately metal-deficient AGB stars should already have a 
substantial DUP parameter of A ~ 0.5. The fact that our models show larger A is generally 
consistent with this trend. It is shown that efficient dredge-up can be found even under the 
assumption of small overshooting efficiency, and that a high numerical resolution is required 
for dredge-up predictions. 
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Table 1. Model sequences used in this paper 



track ID 


/PDCZ 


/CE 


M ZAMS /M 


Z 


E79-D4 a 


0.002 


0.016 


5 


0.0001 


E79-D9 


0.016 


0.016 


5 


0.0001 


E79-D10,ll 


0.002 


> 0.032 


5 


0.0001 


E37-D2 


0.016 


0.016 


5 


0.02 


E85-D2 


0.016 


0.008 


4 


0.0001 



"benchmark sequence 
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Fig. 6 - Panel (1): Evolution of the La- 
grangian location of the H- and He-burning 
shells and the convectively unstable layers (dark 
shade = PDCZ, light shade = envelope con- 
vection). Panel (2-5): Abundance profiles 
at times indicated by thin vertical lines in 
panel(l). Labelled regions of special interest 
are discussed in the text. A blow-up of region 
F in panel (iv) is shown in Fig. 7. 



5? 3.6 




Fig. 1.— HRD of pre-AGB evolution. The 
plus signs along the track indicate selected 
points along an evolution sequence computed 
by Girardi et al. (1996) with the same initial 
mass and metallicity. 



M=5M sl 
Z=0.0001 
sequence E79-D4 




40 60 

t-t / 1 0OOyr 



Fig. 2. — Evolution of the hydrogen- free core 
(solid line, defined as M H = m r (X H < 0.37)) 
and the helium-free core (dashed line, de- 
fined as M He = m r (X He < 0.49)) during the 
TP-AGB. t = 99385139.9 yr at maximum of 
L Hc during first thermal pulse. M = 5M , 
Z=0.0001, sequence E79-D4. See text for an 
explanation of the small DUP at t — to ~ 
15 • 10 3 yr. 
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Fig. 3. — Evolution of the hydrogen-burning 
luminosity (top panel) and the helium burn- 
ing luminosity (bottom panel) for the same 
sequence as shown in Fig. 2. 



Fig. 4. — Evolution of the stellar radius 
(top panel) and the stellar luminosity (bot- 
tom panel) for the same sequence as shown 
in Fig. 2. 
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0.9645 




0.961 



18.12 



18.13 



18.14 



18.15 
t / 1 OOOyr 



18.16 



18.17 



18.18 



Fig. 5. — Evolution of the hydrogen-free 
core, the helium-free core, and the convec- 
tion zones during the sixth TP ( M = 5M Q , 
Z=0.0001, sequence E79-D4). Except for a 
short period immediately after the He-flash 
at t = 18.14 • 10 3 yr, the bottom of the en- 
velope convection zone and the H-free core 
coincide in this representation. Note that the 
third DUP penetrates into the He-free core. 
to = 99385139.9 yr at maximum of L^ e dur- 
ing first thermal pulse. 
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Fig. 6. — for caption see page 25 
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5e-07 
m r -0.959552/M„, 



1e-06 



1 .5e-06 



Fig. 7. — Abundance profiles of species rel- 
evant for the s-process in region F of Fig. 6, 
panel (iv). The region on the left is the top- 
layer of the C/O core, the region on the right 
is the bottom of the envelope convection zone 
at the end of the third DUP after the tenth 
TP. The helium abundance in the layer be- 
low the convection zone is low because the 
third DUP has reached below the He-shell 
(see text). 
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Fig. 8. — Profiles of protons, 12 C, nuclear en- 
ergy generation e (left scale), and tempera- 
ture (right scale) at the bottom of the con- 
vert ive envelope during the DUP of sequence 
E85 (Z = 0.0001). The abscissa has been 
shifted by the mass coordinate of the enve- 
lope convection boundary, and the region to 
the right of the vertical line is the convective 
envelope. 
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Fig. 9. — Envelope abundance evolution of 
relevant CNO isotopes for the 5M bench- 
mark case (E79-D4). Two irregular be- 
haviours can be seen. At t ~ 15 • 10 3 yr, the 
DUP is smaller than after the preceding and 
the following thermal pulse as the result of 
insufficient numerical resolution (see §4.2 for 
details). At t ~ 62 • 10 3 yr, some data files for 
the plot routine were lost, and the HBB phase 
following the thermal pulse - most notably of 
12 C - is not plotted accurately. 
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Fig. 10. — Oxygen and sodium abundances in 
globular cluster stars (Kraft et al. 1997; Sne- 
den et al. 1997; Shetrone 1996; Kraft et al. 
1992) and AGB model predictions (Table 1). 
The 5 M Q sequence E79-D4 has been com- 
puted without mass loss and achieves high 
HBB temperatures, while the 4 M sequence 
E85-D2 has been computed with mass loss 
and neither oxygen nor sodium are destroyed 
by HBB. Both sequences start with low O 
and Na abundances indicated by a filled cir- 
cle (see text for details). For the conversion of 
the model abundance mass fractions into the 
logarithmic abundance ratios with respect to 
solar log 16 O = -2.018, log 23 Na = -4.476 
and log 56 Fe/ 56 Fe = -2.3 has been used. 
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Fig. 11.— Abundance evolution in units Fig. 12.— Abundance evolution in the same 
of [X/Fe]=log(X/Fe)-log(X /Fe ), with units as in Fig. 11 for test calculation E79- 
log 16 O = -2.018, log 12 C = -2.479, D10 after the fifth TP (see text for details). 
log 14 N = -2.957 and log 56 Fe/ 56 Fe = t = yr at maximum L He during the flash. 
—2.3. Top panel: sequence D4 

(/pdcz = 0.002, /ce = 0.016); bottom 
panel: case D9 (/pdcz = /ce = 0.016) (see 
Table 1). 



